NONCOMMUTATIVE BESSEL SYMMETRIC FUNCTIONS 



JEAN-CHRISTOPHE NOVELLI AND JEAN- YVES THIBON 

Abstract. The consideration of tensor products of O-Hecke algebra modules leads 
to natural analogs of the Bessel J-functions in the algebra of noncommutative 
symmetric functions. This provides a simple explanation of various combinatorial 
properties of Bessel functions. 



1. Introduction 

It is known that the theory of noncommutative symmetric functions and quasi- 
symmetric functions is related to O-Hecke algebras in the same way as ordinary sym- 
metric functions are related to symmetric groups. Thus, one may expect that natural 
questions about representations of O-Hecke algebras lead to the introduction of inter- 
esting families of noncommutative symmetric functions. By "interesting", one may 
mean "noncommutative analogs" of the Frobenius characteristics of representations 
of symmetric groups based on combinatorial objects, which may themselve give back 
various identities for the ordinary, exponential, or g-exponential generating functions 
of these objects. This amounts to specialize the complete noncommutative symmetric 
functions S n (A) to h n (X), 1, ^ or — !p, respectively. 

Examples of this situation can be found in [TZj, where the analysis of the rep- 
resentation of H n (0) on parking functions leads naturally to the combinatorics of 
the noncommutative Lagrange inversion formula, and to the introduction of non- 
commutative analogs of various special functions, such as the Abel polynomials, the 
Lambert binomial series or the Eisenstein exponential, and allows one to recover in 
a straightforward and unified way a number of enumerative formulas. 

The present paper addresses the following question. The O-Hecke algebra is the 
algebra of a monoid, hence admits a natural coproduct for which the monoid elements 
are grouplike. This allows one to define the tensor product of O-Hecke modules, which 
induces on quasi- symmetric functions an analog of the internal product of symmetric 
functions. What are the properties of this operation, and of the dual coproduct on 
noncommutative symmetric functions? 

It turns out that the second part of the question is the most interesting. Basically, 
the answer is: the dual coproduct governs the combinatorics of Bessel functions. In- 
deed, its explicitation leads to the introduction of noncommutative analogs J n (^4, B) 
of the J-functions of integer index, of which a few basic properties are readily es- 
tablished. Then, the above mentioned specializations (and other more complicated 
ones) give back various classical enumerative formulas. 
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2. Background 



2.1. Notations. Our notations for noncommutative symmetric functions are as in 
[T2*l IT3j . The Hopf algebra of noncommutative symmetric functions is denoted by 
Sym, or by Sym(A) if we consider the realization in terms of an auxiliary alphabet. 
Bases of Sym n are labelled by compositions I of n. The noncommutative complete 
and elementary functions are denoted by S n and A n , and the notation S 1 means 
Sii • • ■ Si r . The ribbon basis is denoted by Ri. The notation I \= n means that / is a 
composition of n. The conjugate composition is denoted by 7~. 

The graded dual of Sym is QSym (quasi-symmetric functions). The dual basis of 
(S 1 ) is (Mi) (monomial), and that of (Ri) is (Fj). 

The Hecke algebra H n (q) (q G C) is the C-algebra generated by n — 1 elements 
Tx, . . . , T n _i satisfying the braid relations and (Tj — l)(Tj + q) — 0. We are interested 
in the case q = 0, whose representation theory can be described in terms of quasi- 
symmetric functions and noncommutative symmetric functions [To] IHj. 

The Hopf structures on Sym and QSym allows one to extend the A-ring notation 
of ordinary symmetric functions (see ^3], and JHj for background on the original 
commutative version). If A and X totally ordered sets of noncommuting and com- 
muting variables respectively, the noncommutative symmetric functions of XA are 
defined by 



Thanks to the commutative image homomorphism Sym — > Sym, noncommutative 
symmetric functions can be evaluated on any element x of a A-ring, S n (x) being 
S n (x), the n-th symmetric power. Recall that x is said of rank one (resp. binomial) 
if a t (x) = (1 — tx)^ 1 (resp. a t (x) = (1 — t)~ x ). The scalar x — 1 is the only element 
having both properties. We usually consider that our auxiliary variable t is of rank 
one, so that at (A) = a\(tA). 

The argument A of the noncommutative symmetric functions can be a "virtual 
alphabet". This means that, being algebraically independent, the S n can be spe- 
cialized to any sequence a n G A of elements of any associative algebra A. Writing 
ot n = S n (A) defines all the symmetric functions of A, and allows one to use the power- 
ful notations F(nA), F((l — q)A), etc., for more or less complicated transformations 
of the specialized functions. 

The (commutative) specializations A = E, defined by 



(1) 



a t (XA) = Y J t n Sn(XA) = Y[ (Tt*{A) = ^^(X)^) • 



n>0 x£X I 



S n (E) 



1 



and A = — !— , for which 

1 — a 7 



(3) 




are of special importance. 
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2.2. Noncommutative analogs of special functions. Since the discovery by D. 
Andre of a combinatorial interpretation of tangent and secant numbers, several clas- 
sical generating functions have been lifted to the algebra of symmetric functions, 
and more recently, to noncommutative symmetric functions. The general idea is as 
follows. Given the exponential generating function 



of a combinatorial sequence a n e N, one looks for a noncommutative symmetric 
function F(A) such that F(tK) = /(£). The noncommutative analog is interesting 
when F n (A) can be directly interpreted as the formal sum of the combinatorial objects 
counted by c n , under the embedding of Sym into some larger algebra. For example, 
in the case of tangent and secant numbers, the series 



becomes the formal sum of the alternating permutations (shapes (2 ra ) and (2 n l)) 
under the embedding of Sym in FQSym ^2]- One can also find in J2| the noncom- 
mutative Eulerian polynomials, and in ^7], analogs of the Abel polynomials and of 
the Lambert and Eisenstein functions. 

In general, F n turns out to be the characteristic of some projective O-Hecke mod- 
ule. Projective modules are always specializations of generic modules, thus also 
representations of the symmetric group, whose Frobenius characteristic are then the 
commutative images F n (X). In general, setting X = gives back an interesting 
q- analog of f(t). 

In this note, we shall show that the consideration of O-Hecke modules obtained from 
a natural notion of tensor products leads immediately to noncommutative anlogs of 
the Bessel J (or I) functions. Here, we need two alphabets A and B, and we are led 
to the combinaorics of bi-exponential generating functions. 

3. Tensor products of O-Hecke modules 

3.1. The O-Hecke algebra H n (0) is the algebra C[n„], where the monoid H„ is gen- 
erated by elements tti, . . . , 7r n _i (7T» = 1+TA satifying the braid relations 

(6) ITillj = TTjTCi \i — j\ > 1 

(7) -K i 'Ki+\'Ki = 7Tj + i7rj7rj + i 

(8) 

and the idempotency condition 



(4) 




n>0 
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Hence, tensor products of iJ n (0)-modules can be defined, and it is obvious from the 
definition of the simple module Sj that 

(11) S H ® S K = where Des(J) = Des(#) n Des(K) . 

This induces an internal product A on QSym n = G (H n (0)), similar to the internal 
product of symmetric functions, such that 

(12) F H A F K = F 7 where I = H A K, that is, Des(J) = Des(H) n Des(ii~) . 
By duality, this defines a coproduct on Sym n , given by 

(13) 7a#/= Rh®Rk- 

Des(7)=Des(H)nDes(A") 

3.2. There is a canonical involution i on H n (0), defined by 

(14) L(7[) = 7fj = 1 - TTi , 

so that we can regard H n (0) as C[II n ] as well. Hence, we have another tensor product, 
defined from the coproduct 

(15) <5 V 7T = 7f (g) 7f for 7T G n n , 

which induces a second internal product V on QSym, 

(16) F H y F K = F 7 where I = HVK, that is, Des(J) = Des(# ) U Des(fsT) . 

It is of course sufficient to study one of them. However, it is interesting to observe 
that this second product appears in another guise in jH] , in the process of calculating a 
basis of primitive elements of FQSym. Let us recall this construction. Let p n denote 
the projection onto the homogeneous component FQSym n of FQSym, and let fi q : 
F Q cg> F^ i— > F aUlg p[k] be the multiplication map of FQSym g . The g-convolution of 
two graded linear endomorphisms /, g of FQSym is defined by 

(17) / O q g = ^ q o(f®g)o A. 

For q = 1, this reduces to ordinary convolution. We are interested in the case q = 0. 
For a composition I = (z 1; . . . , i m ), let 

(18) Pi = Ph®w-®oPi m - 

It is proved in |H] that the pi are mutually commuting projectors, and more precisely 
that 

if|/|^|J|. 
pi\/j otherwise . 

Hence, j : Fj i— > pj defines an embedding of (QSym, V) in the composition algebra 
of graded endomorphisms of FQSym. Moreover, 

(19) 7T=X)(-1) ,(J) " 1 P/ 

IJI>1 



Pi°Pj 
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which is a projector onto the primitive Lie algebra of FQSym, is the image of the 
primitive element ^2 n M n of QSym under j, and it easy to see that more generally, 
for any / G QSym 

(20) (j ® ]){A QSym f) = A FQSym o j(f) . 

However, j does not map the usual (external) product of QSym to the ordinary 
convolution of endomorphisms. It is nevertheless interesting to pull back the 0- 
convolution to QSym, by defining 

(21) F I Q Fj = F I .j, 

where I ■ J means as usual concatenation of the compositions. Then, we have a 
splitting formula 

(22) (A 0o h 0o • • • ©o fr) vg = MCA - • • fr) v r A r QSym (g)] 

analogous to the one satisfied in Sym. 

It can be shown that the involution i maps the simple module Sj and the inde- 
composable projective module Pj to S/~ and P/~, respectively. 

3.3. Identifying as usual a tensor product F <g> G with F(A)G(B), where A and B 
are two mutually commuting alphabets, we have 

(23) a^XA) A a x (XB) = ^ F k (X) 1a (R k ) = ^(XA) , 

K 

which may be compared with the following identity relating the internal product * 
of Sym and its dual coproduct 5F = F(XY) on QSym: 

(24) o x (XA) * o-i (Y A) = (7i (XYA) = Sax (XA) . 

Theorem 3.1. The coproduct 7 A is a morphism for the ordinary (outer) product of 
non commutative symmetric functions, that is 

(25) lA (FG)= lA (F) lA (G) 

In particular, it is completely determined by the images of the elementary functions, 
7 A A n = A n £g> A n , which implies the combinatorial inversion formula 

(26) (^(-lfA^A^ = Rh®Rk- 

\n>0 J Des(H)nDes(K)=9 

Proof - This is equivalent to Theorem 14.11 below. B 



As we will see, this simple identity has many interesting enumerative corollaries. 
Applying the involution u on the second factor gives the inverse of 

(27) ^(-lfA n ®^ n ") = Yl Rh®Rk~- 

\n>0 / Des(H)nDes(K)=0 

The right hand side of this equality occurs in [T^j, where it is interpreted as the 
decomposition of the algebra H© ra as a bimodule over itself. The inverse of the 
left hand side can legitimately be considered as a noncommutative analog of the 
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Bessel function J , as if we specialize both sides to xE, we recover J (2x). Moreover, 
specializing A to x/(l — q) gives a classical g-analog of J , and the other ones are 
obtained by simple transformations. This first step being granted, it is not difficult 
to guess the correct definition of the noncommutative analogues of the other J v . This 
will be done in the forthcoming section. 

4. Noncommutative Bessel functions 

4.1. Let A and B be two mutually commuting alphabets. The noncommutative 
Bessel functions 3 n (A, B) are defined by the generating series 

(28) zHJ ^ B ) = ^-i/z{A)a z {B), 
that is, 

(29) J n (A, B) = J2(-l) m ^m-n(A)S m (B) . 

m>0 

For A = B = xE, this is the usual Bessel function J n (2x). In particular, 

(30) MAB) = ^(-l) m A m (A)5 m (5) 

m>0 

can be regarded as A_i(J), for the virtual alphabet J = (A, B) such that 

(31) A n (J) = A n (A)S n (B) . 
This defines an embedding of algebras 

j : Sym -> Sym(A, B) = Sym (g) Sym 
A n (A) ^ A n (J) = A n ® S n . 

It is not difficult to describe the image of the ribbon basis under this embedding. 
We need the following piece of notation. For two compositions / and J of the same 
integer n, we define the composition K = I\J of n by the condition 

(33) Des(K) = Des(/)\ Des(J) (set difference) . 
Then, we can state: 

Theorem 4.1. The image of Rk by j is 

(34) R K (S) = J2 RM)Rj(B) . 

I\J=K 

Proof - The formula is true for K = (l n ) by definition. The general case follows by 
induction on l(K~), the number of columns of the ribbon diagram of K. Indeed, it 
suffices to prove that 

(35) Rk{T)Ri™{T) = iW(J) + iW«(J) , 

which follows easily from the usual multiplication rule of ribbon functions. B 
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Corollary 4.2 ([3J. Let a n be defined by 



1 y-n 

n>0 



Then, a n is equal to the number of pairs of permutations (er, r) G & n x © n st/c/i i/iai 
Des(a) C Des(r). 

Let 5 be the linear operator on Sym (acting on the right) defined by 

g(il,-,V) rj_ g(il,...,i r -l) 

It has the following properties (see [Hj, Prop. 9.1): 

(38) (FG) 9= F • (G d) + (F d) ■ G , 
where G denotes the constant term of G, and 

(39) R, 9= {^..-- ''{;:> = \ . 
In particular, if Go = 0, 

(40) (1-G)- 1 d= (l-G)-\Gd). 
Let us apply this with d=ds acting only on Sym(£>) to 

(41) J (A BY 1 =(l- J2(-^ 1 MB)S n (A)\ = S I (A)R I (B) . 

\ n>l / I 

We obtain 

(42) 3 (A,B)- 1 3„ 1 (A,B) = Y,S I {A){R I d (B)) 

i 

Corollary 4.3 ([3J. The coefficient c n in 

(43) MM = Tc ^ 



is equal to the number of pairs of permutations (a, 0) G 6^ such that Des(a) C Des(/3) 
and P(n) = n. 

4.2. Bessel-Carlitz functors. Let F be the functor which associates with a pair of 
vector spaces (V, W) the graded subalgebra of the exterior algebra A(V © W) 

(44) F( V, W) = A n (V) ® A n (W) . 

n>0 
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This is a quadratic algebra (see|18j). If (t>j), (vjj) are bases of V and of W, the 
relations are as follows. For i < k and j < I, 



(45) 


i k 


+ 


ki 
J 1 . 


= 0, 


(46) 


i k 
J 1 . 


+ 


i k 


= 0, 


(47) 






a 

J 1 . 


= 0, 


(48) 






i k 
3 3 


= 0, 



where 



Hence, the Koszul dual G(V, W) 
presented by 



(49) 



i k 




ki 




i k 


Jl 




Jl 







F(V, W)- is the quadratic algebra on V* <g> W 
for i < k and j < I. 



The combinatorial investigation of Bessel functions has been initiated by Carlitz 
0. Hence, the polynomial bi-functors defined by F and G can appropriately be 
called Bessel- Carlitz functors. One or two occurences of A can be replaced by S in 
the definition of F. In the mixed case A Cg> S, the best interpretation is probably as 
functors defined on super (i.e., Z 2 -graded) vector spaces V — Vq © V\. 

5. The ^-specialization 

This section is devoted to the interpretation of a few formulas from [3J El i n 
terms of noncommutative symmetric functions. 

5.1. Let 9 C A x A be any binary relation. We denote by 9 the complement of 9 in 
A x A and set 

X =X(A; 9) = {w = ai ■ ■ ■ a n 6 A*\a 1 9a 2 9 . . . 9a n } , 
( ' ,U) Y =Y(A;9) =X(A;9), 

where we write a9b for (a, b) G 9. Note that the empty word 1 and the letters belong 
to both X and Y. 

The ^-specialization Sym(A; 9) is then defined by specifying the elementary sym- 
metric functions 

(51) A n (A;9)= w - 

wexnA n 

The following basic lemma, implicit in |4., generalizes the case 9 = {(a,b)\a > b}. 

Lemma 5.1 (Carlitz-Koszul duality for alphabets). The complete symmetric func- 
tions S n (A;9) are given by 

(52) S n (A] 9) = A n (A;9) . 
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More generally, if one denotes by 9Ad}(w) = {i\ai9a i+1 } the 6-adjacency set of w = 
a i^2 • ■ ■ &n, and by Cg(w) the associated composition of n, one has 

(53) R I (A;8)= y w - 

C (w)=I 

Proof - We need to prove that 

n 

(54) y(-l) k A k (A, B)K-k{A, 9) = 

fc=0 

for n > 0. Let w = uv be such that u G Ak(A,9) and v G A n _k{A, 9). Then if 
last(u)#first(i>), w appears in Ak+i(A, 9)A n _k-i(A, 9), and similarly, if last(u)^first(v), 
then w appears in Ak-i(A, 9)A n ^k + i(A, 9). Moreover, w cannot appear in any other 
product, so that its coefficient in the sum is 0. B 

5.2. The #-Eulerian polynomials. Recall from [12J that the noncommutative Euler- 
ian polynomials 

(55) A n (t;A) = Y,t l(I) Ri(A) 

It=n 

admit the generating function 

(56) ^^^M)=r^jj 

n>0 1 tV ' 

(see |7] for the commutative version of this identity), and since l(Cg(w)) = 6 l adj(w) + l, 
we have immediately 

(57) y t ead ^ +1 w = - — 1 ~ t . A - . 

i-t«i-t(A;9) 

Note that 9&d}(w) + 9adj(w) = n — 1. Replacing 9 by 9, A by t^A, then t by t -1 , 
and simplifying by (1 — t) the resulting expression, we obtain 

(58) V t 9&d ^W = = 1 ~ —J- — r — - — , 

which is Theorem 2 of jllj . 

For a letter c G A, denote by dc the linear operator defined by 

*~ (u if w = uc for some u , 

( 59 ) Wdc= \0 otherwise. 

Then, as in (|37|). for any series F without constant term, 

(60) (l-F)- 1 dc=(l-F)- 1 -(Fd c )- 
The same is true for the operators 

(61) D c = dc -c 
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where C is a subset of A. Applying this to (JHEJ), we obtain 

( 62 ) E t9adKW)w = i "v W€XC(t F^P • 

«fco 1 -i)' (,o) ^ 

which is Theorem 3 of 

5.3. The #-Major index. If one defines the #-Major index by 

(63) 0maj(u>) = i 

one has clearly 

(64) Yl q 9mai(w) ™ = E g maj(/) iM^; f) = (q) n S n (-^-; e) 

weA n Ihn ^ " / 

where as usual 



(65) 



n>0 

6. Double Eulerian polynomials and Bessel functions 



6.1. The noncommutative Bessel function Jo (A, B) can now be properly interpreted 
as a generating series of 9 elementary symmetric functions, if we interpret J as the 
product alphabet A x B, endowed with the relation 

(66) (a, b)9{a', V) <S> a > a' and b < b' . 

As is customary, we denote words over A x B by biwords 



(67) w = [u, v] 

Observing that 



ii g A n , j/eF 



Des(w) nDes(w) = Des(w)\ Des(i> ) 



(68) flAdj 
we can now write 

\p -0adj (w) l(w) = 1 ~ j 

^ J ((l — t)z; A,B)—t 

( 69 ) w=(u,v)e(AxB)* uu 1 ' ' > 

= J2z lKl t l{K) - 1 R K (A 1 B; 



K 

where from now on we shall use the notation 



(70) J (x; A, B) = A_ X (J) = \„ X {A, B; 9) . 

The coefficient of z n is the nth double Eulerian polynomial, denoted by A n (t; A, B; 9). 
Setting A = B = E, we recover the enumeration of pairs of permutations (a,/3) G 
& n x & n by the cardinality of Des(a) R Des(/?) (cf. [3]). 
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7. The Fedou-Rawlings polynomials 

By considering simultaneously the specializations of (jnHj) to all positive q and p- 
integers, A4 = [i + l] q and Bj = [j + l] p , one arrives at the five parameter gen- 
eralizations of the double Eulerian polynomials introduced by Fedou and Rawlings 
P. 

For w G A n , where A is the infinite chain A = {ai < a-i < . . .}, let q w be the image 
of w by the multiplicative homomorphism a, 1— > c/ _1 . Writing, for a composition 7 of 

n 

(71) i2,(A) = E E ™ 

C(cr)=/ Std(tl>)=<7 

and taking into account the identity 

™des (a ~ 1 ) „coimaj (cr) 

(72) V x max(M,) g u ' = — 

^ (xq;q) n 

where coimaj(<r) denotes the co-major index of cr -1 , 

(73) coimaj (cr) = E (n — d) , 

dGDesO- 1 ) 

(indeed, it is easily checked that the minimal word v for the lexicographic order such 

that Std(v) = o satisfies q v = g coima jW) ; W e find 

(74) 

V x'iUl, 9 , . . . , q 1 ) = — V X max (™)g«< = 1 V" ^desCa-^^coimajCa) 

i>0 C(u>)=7 V ' ^ n+i C(<t)=J 

so that finally, we recover the double generating series of ^1] 
(75) E^V — 



M .> J ((i-t)^;A,^)-t 

_ ^desris(a,/3) 3 ,des(a" 1 )^des(/3~ 1 )^coimaj(a)pCoimaj(/3) 

^ (x;q) n+1 (y;p) n+1 

where desris(a, (3) = | Des(a)\ Des(/?)|. 

The second generating series of pi] is recovered in the same way. If we denote by 
bj the greatest letter of Bj, then, on the one hand, 

(76) S n (Bj) dbj— S n -i(Bj) . 
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On the other hand, 



(77) ^Rj(Bj) = J- J2 

J>° y C{a) = J 



max(v) v 



max(w) = last(«) 

1 

y 



C(cr) = J Std(«)=<r 
cr(n) = n 



max(v)^ 



so that, applying the operator 5^ -fy to the coefficient of x l ]p in (|75|). we obtain 
(78) £ ^-^xV^^)^^')-^-^ 

if I\J=K ij>0 

= ^ rrV f 1 " 5> n ^ - l) n " 1 An(A,,)5 n ( J B i )') ■ S 6 , -b 3 

i,j>0 \ n>l / 

= ^ < j (~ En>l ~ l)- 1 An(^)g n _ 1 (i? J )6 J ) (1 - t) 

3 J_!((l-t)z; Aj, 



X2/ J (^(l - t); A,, ^) - t ' 



Specializing Aj = [z + 1] 9 , 5j = [j + l] p , this becomes, in the notation of jTT], 

m>o J c i ((i- 

_ ^desris(a,/3) a ,des(o- 1 )^des(/3- 1 )^coimaj(a)pCoimaj(^) 

^ (a?;g)„+i(y;p)„ 

n - u a,p G 6„ 

/3(n) = n 

which is equivalent to [TTJ (3)]. Here, 

(80) q,p) := (-l)"J v («[i + l] g , [j + l] p ) . 

The other results of ^T] can be rederived in the same way, by changing the special- 
izations of Ai and Bj. 



8. Heaps of segments and polyominos 

Bessel functions and their multiparameter analogs play a crucial role in the enumer- 
ative theory of polyominos . Elegant combinatorial proofs of such enumerative 
results can be achieved by means of Viennot's theory of heaps of segments [T9~|H]. As 
we shall see, this can also be conveniently formulated in terms of #-noncommutative 
symmetric functions. 
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8.1. A parallelogram (or staircase) polyomino P, which is also the same as a con- 
nected skew Young diagram, can be encoded as a biword 



(81) 



k • 


' 1"n 




U 


J 1 ' 


' Jn 




V 



where j& is the heigth of the /cth column C*., and z^. is the number of common rows 
between Ck and C^+i (with a conventional value i n — 1 for the last column). For 
example, the following polyomino 



is encoded by the biword 



2122111 
2323212 



The biwords corresponding to polyominos are the words over the alphabet 

(83) A = {a ij \i<j} 
satisfying the ^-adjacency conditions 

(84) a ikjk 9 a ik+ljk+1 i k < j k+1 
and the ending condition 

(85) i n = 1 . 

Hence, the generating series (by length) of all biwords satisfying (J%4*|) is 



(86) \ t (A,9) = [\_ t (A, 



i-i 



tl ■ ■ -in 
31 ■ ■ ■ 3n 



i-Et- 1 ) 71-1 * 71 E 

n>l i k >j k+1 •-' 

and restriction of the series to the biwords satifying (|8*K|) is achieved as above by 
applying the operator 



(87) 



D = E 9« 



i>i 



so that we end up once more with a series of the form 
8) 

1 -E(- 1 ) n " lr E 



n>l 









h ■ 




)"( 


J 1 ' 


■■jn 





E(-i) n ^t n E 

«fc>ife+i;«n=i 



n>l 



jl ■■■jn 
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which acquires the structure J\/ Jq once A is specialized to 

(89) = ay'-y , 

the generating series by width, height and area. 

8.2. This can of course be interpreted in terms of heaps of segments. A segment is 
an interval of N*. To each segment, we associate a variable 



(90) a. 



ij 



in our A = {a,ij\i < j}. The monoid of heaps is the quotient of the free monoid A* 
by the commutation relations 

(91) aijOM = a kl aij if j < k 

which means that the segments do not overlap and can be vertically slided indepen- 
dently of each other. 

The first basic lemma of the theory (which is also a special case of the Cartier- 
Foata formula for the Moebius functions of free partially commutative monoids [3]) 
amounts to the calculation of S n (A, 9) for the relation defined by 

(92) dij 6 a M i < I . 

Indeed, with this choice, A n (A, 9) is the formal sum of trivial heaps (products of 
mutually commuting segments arranged in decreasing order), and A n (A,8) is the 
sum of all biwords 



(93) w = a il j 1 ■ ■ ■ a 



?"n3n 



such that < jk+i for all k, those encoding polyominos. We have therefore shown 
that each heap, or, equivalently, each element of A* j = has a unique representative 
of this form. 
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